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We examine entanglement dynamics via concurrence among four two-state systems labeled 
A, a, B, b. The four systems are arranged on an addressable "lattice" in such a way that A 
and a at one location labeled Aa can interact with each other via excitation exchange, and the 
same for B and b at location Bb. The Aa location is prepared entangled with the Bb location, but 
their mutual complete isolation prevents interaction in the interval between actions of an external 
addressing agent. There are six pairwise concurrences on the lattice, and we follow their evolution 
in the interval between external actions. We show how entanglement evolves and may exhibit the 
non-analytic effect termed entanglement sudden death (ESD), with periodic recovery. These loss 
and gain processes may be interpreted as entanglement transfer between the subsystems. 
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I. INTRODUCTION 

One of the challenges of realizing quantum informa- 
tion processing is control of evolution of qubits in the 
presence of environmental noises and manipulation in- 
accuracies [l|- This control may be easily achieved for 
a single qubit but for many qubits, entanglement 
dynamics has been a difficult subject and has attracted 
extensive interest recently rang ing from two-qubit sys- 
tems S II II i, 0, II llj3jllL to continuous variables 



[121 1131 . spin systems 14 L 3" Tal and multi-partite sys - 
tems|l3, El, El [H M, [H M, El . Moreover, in [H, H|, 
proposals have been made for the direct measurement of 
finite-time disentanglement in cavity QED, and real-time 
detection of entanglement sudden death (ESD) has been 
reported very recently [27| . 

The behavior of qubits in quantum computing envi- 
ronments is expected to be determined by external oper- 
ations, addressed to specific qubit locations, that manip- 
ulate transitions (operate gates) in sequences determined 
by quantum algorithms. Within a qubit network at dif- 
ferent locations there may be lengthy intervals between 
such externally mandated gate operations, and the small 
"network" we examine here is intended to provide insight 
as to the preservation or degradation of entanglement be- 
tween remote qubits on the lattice in these intervals. Our 
network has a dual atom-photon basis in the sense that 
each active lattice site incorporates a photonic and a non- 
photonic element that communicate with each other by 
an emission-absorption interaction. A cavity QED inter- 
pretation is possible [3, [28| . 

One should not expect changes in entanglement be- 
tween two arbitrary units to be determined only by their 
mutual interaction strength. Results obtained below 
demonstrate changes in the entanglement of remote units 



of the lattice that have zero mutual interaction and no 
communication, quantum or classical. 

We begin by specifying the interactions that will be 
allowed. These are intended to permit systems A and a 
to exchange a degree of excitation (a spin flip or a pho- 
ton, for example). The same should be true of B and b. 
Thinking in the cavity QED context, a suitable begin- 
ning would be a sum of two separate Jaynes-Cummings 
Hamiltonians (with H = 1) [HIHHO]: 



A 



+ 5(a t crf + cr^a) + wo'a 



^of + g(bW+<r»b) + u>tfb, (1) 

where luq and ui are the frequencies of atoms and cavities, 
respectively, and g is the coupling constant between the 
atoms and their cavities. Physical realizations of our sce- 
nario does not appear out of the question, as the Jaynes- 
Cummings model has been realized in the laboratory in 
several well-known ways 

Clearly the remoteness of site Aa from site Bb is per- 
fect in this case, since there can be no interaction between 
the sites at all with the Hamiltonian above, although a 
familiar type of interaction will take place on each site. 
It is not difficult to imagine the types of external probes 
that could be coupled to the sites, but we will be in- 
terested in evolution of entanglement in a time interval 
between such external controls. 

The eigenstates of this Hamiltonian are products of the 
eigenstates (dressed states) of the separate JC systems, 
which are well known [29| . We will write for either Aa or 
Bb 



(2) 
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We will think of the non-photonic (atomic or spin) states 
as "excited" and "ground" states denoted by |e) and \g), 
and denote the cavity mode's photon states by the pho- 
ton number n. Then the single-atom eigenvalues are 
given by: 

nc + i(A± ^A 2 + G2), (3) 



A, 
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and the dressed states are these: 

lV>o) = \g;Q) (4) 

\4>+) = c n \e;n~ 1) + s n \g;n) (n > 0) (5) 
\i>n) = -s n \e;n - 1) + \c n \g;n) (n > 0). (6) 

In these equations we have introduced some convenient 
abbreviations: 



cos(<9„/2) and s n = sin(0„/2), 



(7) 



where the rotation angle 9 n can be identified with the 
Bloch sphere polar angle and is defined in the usual way: 



A 



cos( 



and sin#„ = 



Gn 



where A = lj — luq is the detuning and 
G n = 2gy/n. 



(8) 



(9) 



In order to achieve a four-qubit model, we restrict the 
cavities to a flip-flip operation. In the examples to be 
treated we will need the true ground state |-0o) = \9\ 0) 
and the two dressed states for n — 1. These three states 
are closed under the JC Hamiltonian for each site. In 
other words, we use only n — 1 in the equations above, 
so the subscript n can mostly be ignored and we will 
frequently drop it (A n — > A, c n — ► c, etc.). 

The paper is organized as follows: After a brief in- 
troduction to Concurrence chosen as the measure of en- 
tanglement in Sec. HH the evolution of two important 
sets of entangled pure states toward decoherence is dis- 
cussed In Sec. IIIII We examine evolution of concurrence 
between all qubit pairs in this four-qubit model. In par- 
ticular, we show the existence of parameter zones where 
entanglement sudden death (ESD) takes place periodi- 
cally. Moreover, we demonstrate a case where the loss of 
atom-atom entanglement is compensated by concurrence 
gain in the photon-photon space. This loss and gain may 
be interpreted as entanglement transfer between atomic 
and photonic variables. The interconnections between 
pairwise concurrences in this four-partite system are dis- 
cussed in detail. We conclude in Sec. HVl 



II. MEASURE OF ENTANGLEMENT 

The JC dressed states are already atom-photon entan- 
gled states, and this entanglement has interesting conse- 
quences if the cavities are fairly highly excited, as Gea- 
Banacloche originally pointed out [34( , but this is a one- 
site entanglement arising simply from the allowed inter- 
action. Here we are only interested in entanglement as a 
non-local phenomenon that bridges any remote separa- 
tion of sites on the lattice. 

In the general context of entanglement we note that 
there is no accepted and practically workable criterion 
for measuring entanglement (for determining separabil- 
ity) of arbitrary four-particle states. Our purposes will 



be satisfied by working with two-particle mixed states 
obtained from the time-evolving four-particle pure state, 
since the two-qubit domain of entanglement is where all 
familiar measures agree. That is, entropy of formation, 
Schmidt number, tangle, negativity, and concurrence are 
numerically somewhat different, but in the two-qubit do- 
mains of their applicability they are in full agreement 
when they signal complete separability, i.e., the lack of 
entanglement. 

We will adopt Wootters' concurrence [35[ as our mea- 
sure in this discussion, mainly for its relevance for mixed 
states and the convenience of its definition and normal- 
ization: 



CO) = max{0, VAi - V X 2 - \A 3 - V A 4 }, (10) 

where the quantities Ai are the eigenvalues in decreasing 
order of the matrix 



C = p((Jy ® a y )p*(a y ® a y ), 



(11) 



where p* denotes the complex conjugation of p in the 
standard basis and a y is the Pauli matrix expressed in 
the same basis as: 



-i 

1 



(12) 



Notice that 1 > C > 0, where C — indicates zero 
entanglement and C — 1 means maximal pure state en- 
tanglement, as in a Bell state, for example. 

It will turn out that several universal features will ap- 
pear in the entanglements arising from our four-qubit 
model. The simplest ones are these. First, all reductions 
to two-qubit form, obtained by tracing over the other two 
qubits, will yield a two-qubit mixed state always having 
the "X" form [H: 



P = 



a z 

6 

c 

z* d 



(13) 



where a + b + c + d= 1. Second, since the concurrence 
of this mixed state is easily found to be 

C = 2max{0, |z|-v / ad} = 2max{0,Q}, (14) 

it is clear that Q, defined as 

Q = \z\-Vad, (15) 

will be an important quantity. We will mention at the 
end certain "conservation" properties that derive from 
Q in some cases because it can be negative, whereas C 
cannot. 



III. PAIRWISE ENTANGLEMENT DYNAMICS 
OF QUBITS 

The entanglement dynamics of qubits is generally very 
complex even for the case of two qubits, as decoherence 



3 



inevitably evolves a pure state into a mixed state. It is 
easy to see that for an arbitrary initial state of the two- 
atom system there is no exact solution to allow a clear-cut 
analysis. However, we find that two sets of pure initial 
states for two atoms A, B can be examined exactly in 
our model. These initial states are superpositions of the 
two types of Bell states: |$ ± - 1 ) ~ \eA,e B ) ± \9a,9b) and 
\^f^) ~ le^, <7b) ± |<7A,es) ■ Thus we denote superposi- 
tions within each type as follows: 



\®Ab) = coaa\eA,e B ) + sma\gA,9B) 



and 



\^A1 



cosa\e Al g B ) + sma\g A ,e B ) 



(16) 



(17) 



These initial states are obviously interesting and physi- 
cally relevant. It is easy to see that a — ir/i reproduces 
the standard Bell states. It should be noted for certain 
types of initial mixed states such as X-states (3(| , explicit 
results can also be obtained. 

In the two-site situation under consideration there are, 
in principle, six different concurrences that can provide 
information about the bipartite entanglements that may 
arise. With an obvious notation we can denote these 
as C AB , C ab , C Aa , C Bb , C Ab , C Ba . All except C Aa 
and C Bb report remote entanglements across the (Aa)- 
(Bb) gap, and we will be dealing with these four cases. 
That is, in what follows, we will examine in detail the 
entanglement evolution between two atoms, two cavities 
or between one atom and the opposite cavity. 



A. Partially entangled Bell states |3>ab) 

In this section, we will take a combination of the 
Bell states as our initial state for each of two atoms la- 
beled A and B, located without communication in two 
separate cavities with modes labeled a and b. The chosen 
total initial state for the atoms plus cavities is: 

|*(0)) = \<S> AB )®\0 a ,0 b ) (18) 
= (cosa|eA,e B ) + sina\g A , g B )) ® |0 Q ,0 6 ). 

To prepare for the time evolution we express these 
states in terms of the dressed eigenstates given below 
(|4]), and obtain: 

|e Al a ) = c\^t)-s\^) 
IffA, lo> = s\^)+c\ip^) and 
\9A,0 a ) = \tpo)- 



(19) 



Thus the initial atom-atom entangled state is 

|$(0)) = cosa|eA,0 Q ) ® \e B ,0 b ) 

+ sma\g A ,O a ) ® |5b,Q&) 

= cos a{c\^) A - s\^) A ) ® (c\tpt) 



Evolution in time is now easily arranged because the 
evolution of the dressed states is immediate: 



|V ± (*))=e 



|^(0)). 



(21) 



Since the combination of coefficients in our 1^(0)) 
uniquely associates c with \ip + ) and s with \tp~), the time 
evolution can be transferred to the c and s symbols. We 
will henceforth consider them carrying the time-evolution 
exponents. We will use the notation cq and sq to refer to 
their values at t — (no relation to the n = subscripts 
in Eq. ([7|). Then we can write (temporarily indicating 
explicit time dependences for the c's): 

|$(t)) = cos a(c(t)\ip+) A -s(t)\^) A 

[<t)\i>i) B -s(t)\^) B 



(« 



(22) 



where the |V ,± ) states will continue to refer to the states 
at t = 0. 

In order to take traces over individual atoms or cavities 
we need to revert to the bare bases and this leads to: 



= cos a\c(t)(co\e A , a ) + s \g a Aa)) 

- s(t)(-s \eA,0 a + ca\gA,l a ))) ® (c(t)(c o \e B ,0b) 

+ s Q \g B , lb)) - s(t)(-s Q \e B ,0 b + c \g B , l b ))j 

+ sma\g A ,0 a ) ® \g B ,0 b ) 

= cos a(^(cc + ss )\e A ,0 a ) + (cs - sco)\gA, l a )J 

® f(cco + ss o )\e B ,0 b ) + (cs Q - sco)\g B , l b )^j 

+ sina\g A ,0 a )®\g B ,0 b ). (23) 



For orientation we will first focus on just C AB , which 
has been noted previously [28|. This requires that we 
trace out a and b photon modes, in order to get the two- 
qubit mixed state needed for calculation of AB concur- 
rence. The projections that are needed are: 

(0 a ,0&|$(t)) = cos a(cc + ss Q ) 2 \e a ,e B ) +sma\g A ,g B ) 

(l o ,0&|$(t)) = cos a(cs - sco)(cc + ss )\g a, e B ) 

(0 a ,l b \$(t)) = cos a(cc + ssq){cs q - sc )\e A , g B ) 

{laMHt)) = 0. (24) 

These are simple enough to see that the AB mixed 
state has the form mentioned in the previous Section: 



P AB = 



a z 

6 

c 

z* d 



(25) 



-sIV'i 



(20) 



for which the concurrence has the stated form 

C AB =2max{0,\z\-Vte}. (26) 
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We easily find the following 

\z\ 2 = sin 2 a cos 2 a|(cco + sso)| 2 

•2 2/4,4,022 r.\2 

= sm a cos a{c Q + s + 2c s cos St) , 

\z\ — sin a cos (x(cq + s + 2cqSq cos 5t), 

b = c = cos 2 a\cco + sso\ 2 \cso — sco\ 2 

= cos 2 a(cg + Sq + 2cqSq cos <$t) 

x CqSq (2 - 2 cos St). (27) 

For simplicity we will evaluate this in the resonance 
case, 8 n = ir/2, where cq = sq = l/v2- Then we find 

\z\-Vbc= icos 2 a(2+2cos<Si)[tana-sin 2 ((5t/2)], (28) 

where 5|a=o = G, so the expression for concurrence turns 
out to be: 

C AB = 2max{0,Q AS }, (29) 
where Q AB = cos 2 acos 2 (Gi/2)[tana - sin 2 (Gt/2)]. 




FIG. 1: Surface plot of the concurrence C AB for the |*ab) type 
of initial state, as a function of time and the parameter a, from 
Eq. (H) . 

Figs. [TJ O and [3] show that C AB has a novel feature. 
The entanglement non-smoothly becomes and stays zero 
for a finite interval of time for a continuous range of a 
values. This effect was first noted in [a and has been 
found in several contexts ( See, e.g., [j&llH Il3l. I20L |37|). 
It has been referred to as "entanglement sudden death" 
(ESD). Without commenting on the physical implica- 
tions here, this is due simply to the fact that Q AB can 
take negative values. For a < ir/4 there is always a 
time interval in which entanglement remains zero while 
for 7r/4 < a < 3ir/4 entanglement becomes zero only 
momentarily (exactly at multiples of t = tt/G). 

From Eq. (|19[) it can be shown that for resonance, after 
a time interval of n/G, the state |e^,O a ) turns into the 
state \gA, l a ) U P to an overall phase of n/2 and vice versa 
while the state \gA, O a ) does not change. The same is 
true for the B and b parts, of course. Using this fact it 




Gt/n 

FIG. 2: Contour plot of the concurrence C AB for the \<S>ab) type of 
initial state as a function of time and parameter a. The innermost 
contours indicate a concurrence value of 0.9, and the value drops 
by A c = 0.1 between two consecutive contours. Regions of sudden 
death are painted in gray. 




12 3 

Gt/it 

FIG. 3: Plots of atom-atom and cavity-cavity concurrences for 
three values of the superposition parameter a. Entanglement sud- 
den death occurs for all superpositions except the pure Bell states 

with 01= (1± 



can be shown that the |$(t)) of Eq. (|25|) is the same as 
|$(t + tt/G)) under an interchange between atomic and 
photonic states (e^ «-► l a ,9A <-> a ,e s «-* lb, 9b ^> Ob) 
up to a phase ir/2 for resonance. This indicates that 
the state seen from the perspective of atomic part of the 
system at time t is the same as the state that is seen 
from the perspective of photonic part at time t +7t/G so 
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the concurrence between the cavities, C ab , is the same 
as C AB up to a phase difference of tt/G. For the case of 
resonance, we have 

\z\ -Vbc = -cos 2 a(2-2cosGt)[tana-cos 2 (GV2)], 

(30) 

so the expression for concurrence turns out to be: 

C ab = 2max{0,Q ab }, (31) 

where Q ab = cos 2 asin 2 (Gi/2)[tana-cos 2 (Gt/2)]. Thus 
a plot of C ab as a function of a and t would look exactly 
like Fig. Q] except for a phase shift by ir to the left on the 
time axis (see Fig. [3]). 

Now we move to the concurrence calculations between 
opposite atoms and cavities. That is, we will calculate 
C Ab and C Ba and C Aa . 

We start with C Ab . For this, we trace out B and a and 
find, 

(.9s,0 Q $(t)> = cos a(cs - sc )(cc + ss )\e A , l b ) 

+ sin a\g A ,0 b ) 
(e B ,0 o |$(t)) = cos a(cc + ss ) 2 \e A ,O b ) 
(5b>1o|$(*)) = cos a(cs - sc ) 2 |5a, lb) 
(e B ,l 1 <&(*)) = cosa(cs - sc )(cc + ss )|5a, 6 ). 

(32) 

|z| = | sin a cos a | Co so V2 — 2 cos St 

x \Jcq + Sq + 2cqS 2 cos ft 

6 = cos 2 a (cq + Sq + 2c 2 Sq cos ft) 2 

c = cos 2 a c^sl (2- 2 cos ft) 2 . (33) 

For resonance 

Q Ab = |z| -Vte 

= icos 2 a|sinGt| (2| tana| - | sinGiQ. (34) 

Then the concurrence is 

G Ab = 2max{0,Q Ab }. (35) 

Now if we look at the initial state in Eq. (fT5| we see 
that under the transformation A <-> B,a «-> b the state 
remains unchanged. That means C Ba = G" 46 at all times. 
Eq. (f3"4")) suggests that the maximum value C Ab can take 
is less than 1 so the parts A and 6 cannot be maximally 
entangled unlike parts A and B (or a and b). 

As seen from plots (U]) and (O, there exist some pa- 
rameter zones, where entanglement sudden death occurs 
for C Ab in a periodic manner. This interesting feature 
is not expected to be a generic feature for an arbitrary 
initial state, as shown in the next section. 

Similarly, we can find the explicit expressions for C Aa 
and C Bb . The two cavities are not distinguishable, so 




FIG. 4: Surface plot of the concurrence C Ab for the |<3> ab) type 
of initial state as a function of t and a. 



| 0.5 




















0.5 1 1.5 2 

Gt/JT 



FIG. 5: Contour plot of the concurrence C Ab for the \Phi AB ) 
type of initial state as a function of time and parameter a. The 
innermost contours indicate a concurrence value of 0.3, and A c = 
0.05 between contour lines. Regions of sudden death are painted 
in gray. 



we expect these G's to be the same and they are. For 
resonance they become: 

C Aa = |sinGi|cos 2 a = C Bb . (36) 

Before ending this section, we want to stress that the 
onset of entanglement sudden death (entanglement de- 
cays to zero in a finite time and remains zero for at least a 
finite time) is not due to the special mathematical struc- 
ture of Wootters' concurrence. As we have already men- 
tioned, ESD simply signals the finite-time onset of sepa- 
rability (lack of entanglement), and this is a feature that 
will be confirmed by all entanglement measures. More 
explicitly, if a state at time td is separable as measured 
by concurrence, then at the same time it is a separable 
state by all other valid measures - entanglement of for- 
mation, partial transpose negativity, etc. 
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B. Partially entangled Bell States \$>ab) 

Now we examine entanglement dynamics for the fol- 
lowing initial state: 

|*(0)) = \^ab)® |0„,0 6 ) (37) 
= (cosa\e A ,g B ) + sma\g A ,e B )) <8 |0 Q ,0 b ). 



Following the methods of the previous subsection llll A[ 
we find: 

|*(t)) = co8a(c(t)\il>t)A- 8{t)\tl>i) A ) ® K)b 

+ sma^V® (c(t)|^) B -s(t)|^r)s)( 38 ) 
Changing into the bare basis: 

= cosa(c(t)(c |e A ,O a ) + s Q \g A ,l a }) 

- s(t)(-s \e A ,O a ) + c Q \g A ,l a ))j ® |.9s,0 s ) 
+ sina|s-yi,Oa) (8 ^c(t)(c |e B , b ) + s Q \g B , h)) 

- s(t)(-s \e B ,O b ) + Co\g B ,l b ))j. 

(39) 

Again, we start with C AB . The projections that are 
needed are: 

(0 O) 6 |*(*)) = cos a(cc + ss )\e A ,g B ) 

+ sin a(cc + ss )\g A ,e B ) 

(l a ,O b \^(t)) = cos a(cs - sc )\g A ,g B ) 

(0 o , = sin a(cs - sc )\g A , : g B ) 

(l OJ li|*(*)) = 0. (40) 

The density matrix that these projections give is in 
the form of Eq. (j2"5|) again. Yet, this time we have only 
three non-zero diagonal entries instead of the four of the 
previous section. This leads to a zero \fb~c and makes 
it impossible to find sudden death since in this case the 
function Q AB cannot take negative values. 

Continuing the calculations, 

\z\ 2 = sin 2 acos 2 a|(cc + ss )\ 4 

= sin 2 a cos 2 a{c^ + Sq + 2cqSq cos 5t) 2 , 
\z\ = | sina cosa|(co + Sq + 2cqSq cos^t), 



b = 



cos 2 cx\csq — scq\ 2 + sin 2 a\cso — scq\ 2 
\cs - sc Q 1 2 
clsl(2 -2 cos St), 
0. 



For resonance we find 
1 



Q AB 
where 



(41) 



he = j | sin a cos a\ (2 + 2 cos St), (42) 



so the expression for concurrence turns out to be: 

C AB — 2max{0,Q AB }, (43) 

where Q AB = | sinacosaj cos 2 (Gi/2). Since Q > 0, we 
get, 

C AB = 2|sinacosa|cos 2 (Gt/2) 

= |sin2a|cos 2 (G£/2). (44) 
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FIG. 6: Surface plot of the concurrence C AB for the \<t> AB ) type 
of initial state as a function of time and parameter a. 
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FIG. 7: Contour plot of the concurrence C AB for the second type 
of initial state as a function of time and parameter a. The in- 
nermost conturs indicate a concurrence value of 0.9, and the value 
drops by A c = 0.1 between two consecutive contours. Regions of 
sudden death are painted in gray. 

From Figs. © and J7J), as well as ©, we see that C AB 
(also C ab ) lacks the feature of sudden death. The plot 
Fig. [U touches the x— axis only at multiples of t = tt/G. 

It can be shown that Eq. (|3"9"|) has the same symmetry 
as Eq. (|23[) under transformations t — > t + n/G,e A <-> 
K,9A «-> Q ,e B <-> h and g B «-> b , so C ab (t + tt/G) = 
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C AB (t) for resonance. By considering symmetry between 
the atom and cavity, again we have 

C ab = 2Q 

= |sin2a|sin 2 (Gt/2). (45) 




Gt/jt 

FIG. 8: Plots of atom-atom and cavity-cavity concurrences for 
three values of superposition parameter a. Entanglement sudden 
death is now absent for all superpositions although C = is reached 
at a single time for the pure Bell states with a = (1 ± ^)tt/2. 

The remaining concurrence between the atoms and 
cavities are C Ab and C Ba . If we trace out B and a, 

(gB,0a\f(t)) = cos a(cc + ss )\e A ,O b ) 

+ sina(cs - sc )\g A , l b ), 

(e B ,O a |*(t)) = sin a(cc + ss )|3a,0 6 ), 

(g B ,la\^(t)) = cos a(cs - sc )\g a, Ob), 

(e B ,la|*(*)) = 0- ( 46 ) 

Again we only have three nonzero projections which 
means either b or c is zero. In either case, calculating 
only the |io| is sufficient to find out the concurrence. 
For the case of resonance, 

Izl = I sin a cos a | - V2 + 2 cos GV2 - 2 cos Gt. (47) 

Then the concurrence is 

C Ab = 2\z\ 

= 2| siaa cos a|| sin Gt/2||cosGt/2|, 

= | sin a cos a \ \ sin Gt | . (48) 



Looking at Eq. (|37| we see that under a transforma- 
tion A <-> B,a <-* b and cos a <-> sin a the state remains 
unchanged. That means if we just applied the transfor- 
mation cos a <-» sin a to Eq. (|4"B")) we would find C Ba , but 

such a transformation does not change this equation, so 
C Ba = C Ab at all times _ 

Eq. suggests that the maximum value C Ab can 
take is 0.5 so the parts A and b cannot be maximally 
entangled unlike parts A and B (or a and b). 

In the exactly same way, for resonance we get 

C Aa = cos 2 a\ sin Gt\. (49) 

and 

C Bb = sin 2 a | sin Gt |, (50) 

respectively. 




o o 



FIG. 9: Surface plot of the concurrence C Ab for the \^! ab) type 
of initial state as a function of time and parameter a. 

Similar to C AB for the initial state f3T|). the plots © 
and (fit)]) clearly show that there is no sudden death fea- 
ture for C Ab . From the calculations it is evident that 
C AB (t) and C ab (t) have only a phase differ ence of Gt/ir 
between them. This can be explained with excitation 
exchange between the the systems. That is for instance 
from the Hamiltonian one can tell that what is lost from 
A in terms of excitation must go to a so that after half 
a Rabi cycle the two systems (A-B and a-b) exchange 
their initial states. It is also interesting that for the sec- 
ond Bell state the summation of C AB and C ab equals the 
initial concurrence whereas for the first type of Bell state 
it is less than the initial one. 

For the |<I>ab) type of initial state we obtain an ex- 
tra entry in the density matrix (\e A , &b){^Ai &b \ if we are 
looking at the entanglement in the atomic part) that was 
not there for the I^ab) type. Existence of this leads to 
the negative term in the Q function and thus to the sud- 
den death. Without this negative term we would expect 
to get the initial entanglement when we added up the 
entanglements of the atomic and photonic parts of the 
system. However this negative term causes a leak in the 
entanglement. 
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Gt/Ji 



FIG. 10: Contour plot of the concurrence C Ab for the \^ab) 
type of initial state as a function of time and parameter a. The 
innermost loops has a concurrence of 0.4 on them. Concurrence 
decreases linearly as we move to the outer loops with A c = 0.1. No 
region of sudden death is seen. 

IV. CONCLUDING COMMENTS 

We have examined entanglement dynamics of several 
subsystems in a four-qubit model in which two atoms are 
locally coupled to the modes of their cavities. A distinc- 
tive feature in this model is that there are no couplings 
between the two atoms or the two cavities. That is, evo- 
lution of entanglement is a pure information-exchange 
evolution, not a more conventional decoherence process 
involving energy interaction and exchange. 

We have studied the onset of entanglement sudden 
death for two sets of pure initial states. We demonstrate 
that there exist finite-time intervals in which the entan- 
glement measured by concurrence remains zero. Due to 
the absence of interactive decoherence in the pure JC 
model, the lost entanglement can bounce back in a finite 
time. For the l^^s) type of initial state, ESD does not 
occur. Although entanglement decays to zero in a finite 
time, it never remains zero. The occurrence of ESD for 
the l^yis) state is not contradictory, and it is easy to 
see physical distinctions that make the two types of ini- 
tial state nonsymmetric. In the I^ab) initial case two 
photons may later be present at the same time, whereas 
there can never be more than one photon present in the 
\$>ab) case. This said, it is still unclear exactly what 
is the governing distinction. This remains an interesting 
open question in entanglement dynamics. In addition, 
it should be pointed that our main results can be ex- 
tended to the case that the cavity fields are not exactly 
in resonance with the atoms. The qualitative features of 
our conclusions will remain the same, and the detailed 



dynamics in non-resonant situations will simply become 
more complicated. 

In conclusion, we can make three comments: Firstly, 
it is not surprising that entanglement can be generated 
by couplings between atoms and photons. What ap- 
pears here is different. The local atom-cavity couplings 
not only cause the creation of entanglement between an 
atom and its cavity, but also generate non-local entan- 
glement between the cavities if the atom-atom subsystem 
is initially prepared in an entangled state. Moreover, as 
shown explicitly, with aid of entanglement between two 
atoms one may generate non-local entanglement between 
one atom and a remote cavity mode. Secondly, we have 
demonstrated in the case of the |\1/as) type of initial 
state, the loss of entanglement C AB in the atom-atom 
system is instantly compensated by concurrence gain in 
C ab for the cavity-cavity system. In fact, it is easy to 
check that we have 

c ab + C ab = | sin2Q ,| = CQnst ( 51 ) 

We may interpret the entanglement gain in the cavity- 
pair as a process of entanglement transfer from atomic 
variables to photonic variables. However, it may also sug- 
gest some kind of "entanglement conservation". This is 
an open issue that is largely unexplored, and we can't ex- 
pect conservation of entanglement in the sense of dynam- 
ical conservation laws, since entanglement is not defined 
as an observable or represented by a Hermitian operator. 
Interestingly, we can follow the negative values present 
in the concurrence with the aid of the parameter Q AB . 
The following equation holds for all the initial states (fl8|) 
and 1|37[) considered in this paper: 

Q AB + Q ab + 2Q Aa \ tana| — 2Q Ab = |sin2a|, (52) 

where the right hand side is nothing but the initial con- 
currence. This is an issue to be explored elsewhere. Fi- 
nally, our model may suggest an interesting link between 
a system's temporal memory and entanglement genera- 
tion. In the case of the JC model, the memory effect is 
due to the fact that the photons absorbed by the cavity 
modes will be coupled back to the atom in a finite time. 
In a more general context, our research suggests that a 
non-Markovian property may play a crucial role in recov- 
ering entanglement in a quantum dynamical system. 
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